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Abstract. We describe an experiment to measure the mass of the Milky 
Way galaxy. The experiment is based on calculated light travel times along 
orthogonal directions in the Schwarzschild metric of the Galactic center. We 
show that the difference is proportional to the Galactic mass. We apply the 
result to light travel times in a 10cm Michelson type interferometer located on 
Earth. The mass of the Galactic center is shown to contribute 10 — 6 to the flat 
space component of the metric. An experiment is proposed to measure the 
effect. 
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1. Introduction 

What is the mass of the Galaxy? We answer this question by computing, in 
a model metric of the Galaxy, the travel time of light beams in a Michelson type 
interferometer located on Earth's equator. We show that the light travel time is 
proportional to the mass of the Galaxy. The result suggests that an experiment 
can measure the mass of the Galaxy. 

We start with a metric in Boyer-Lindquist coordinates, representing on the equa- 
torial plane, the axially symmetric field of a rotating matter distribution with the 
same symmetry 

(1.1) ds 2 = gttdt 2 + 2g t <t,dtd4> + g rr dr 2 + g de d9 2 + g^dcf) 2 

We will calculate and compare the light travel times for beams traveling along 
the radial, polar and azimuthal axes. The interferometer is small compared to the 
scale factors of the field. 

1.1. Light beam traveling along the <f>— direction. In this section we calculate 
the travel time of a beam of light traveling along the 4>— direction. The beam is 
located on Earth's equator (9 = tt/2), a coordinate distance R from the field center. 
The interferometer orbits about the field center with a constant coordinate speed 
v = RCt, where fl is the coordinate angular speed. The trajectory of light is a null 
geodesic, furthermore, dr = d9 = 0: 

(1.2) ds 2 = = g u dt 2 + 2g t4> dtd4> + g^dcf) 
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We neglect the difference between a null geodesic and an arc of an equatorial 
circle (valid when the coordinate arm length of the interferometer I <C R). The 
world line of the light beam between the two mirrors is: 



-QtA±- >" 2 



J9t4, ~ 9tt94>4> 

(i.3) t± = * {<p-h) 

gtt 

We choose </> = 0. The positive sign is used when the light beam is traveling 
in the same sense as the interferometer. The world line of the end mirror of the 
equatorial arm is: 

(M) 

where $ is the angle subtended by the interferometer at the field center. The light 
ray reaches the mirror when 



-9t<t> + J g 2 ^ - gttg^ 

(L5) Vu * = — 

Solving for the azimuth gives 



(1.6) i 



g t4> n + g tt - SlyJ 5 t 2 - g tt g^ 
and the coordinate time 



5*0 ~ \l9u ~ 9ttg<f>4> 

(1.7) U = 



st^^ + s« - - gttgu 

The world line of the reflected ray is 



-gt<t> - \/ 9td> ~ gttg<t><t> 

(1.8) t = t r + 1 (4>-4>r) 

gtt 

and of the beam splitter is 



The ray returns to the beam splitter when 



-gt<p - \l9tA ~ 9ttg<t><t> 

(1-10) 7i =t r + ~ ^-^) 

Si <7it 

or when the azimuth is 
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gttU + gt<p<j> r + \ 9ts ~ 9tt94>4>4> r 

(1.11) fa = n 

and coordinate time 

(1.12) ti 

9tt + gt<f>^ + Qj9t4, ~ 9ttg<p<f> 

Substituting the expressions for fa and t r gives 



9tt +9t4>^ + fi y 9% - 9tt9U 
gttU + g^fa. + \J 9t4, ~ 9~tt9H>4> r 



U = 2$- 



9t<p ~ 9ttg<p<t> 



< ~ 1 ' 13 ' ) " gtt + Zgt^ + Wg^ 

In terms of proper time of the interferometer this is equivalent to 



(1.14) 
(1.15) 



-\Jgtt + 2gt<t>n + Wg<t,<t>ti 



2 $ / gU_ - gttg^ 

c V g tt + 2g t4> Q + VPg^ 



1.2. Light travel time along the r— direction. We repeat the calculation for a 
light beam traveling inward from the beam-splitter at r = R to the end mirror a 
coordinate distance I at r = R' , where R > R' . In this case d9 = 0. The coordinate 
length of the interferometer arm is I. The world line of the mirror is 



(1.16) 
(1.17) 

The world line of the light ray is now 



r = R' 



R' 



(1.18) 




gtt g 



gu 



> dr 



The negative sign of the radical is used because the light beam is traveling inwards. 
Assume for the moment that ^ is small enough to be kept to first order only. The 
expression simplifies to: 



R' 



(1.19) 



t ~ 



9t<p d<t> 
gtt dr 



gtt 



dr- 



ill the example under consideration we can approximate the space trajectory of 
the light beam as a straight line starting from the origin at (R, 0,0) and passing 
through (R ; , 0, fa). We can relate r with 4> 



4 



MUNAWAR KARIM, ANGELO TARTAGLIA, AND ASHFAQUE H. BOKHARI 



(1-20) r= - , 

cos <p + tan £, sin < 

£ is fixed by imposing 



(1.21) R' 



cos 4> i + tan £ sin <\> i 
This gives 

„„, R-R' cos fa R-R' , 

(1.22) tan £ = — — - y ' w for ^ < 1 

R'smfa R'fa 



(1.23) 



1? 



/ , R—R' COS r/>- 

cos H p , . sin < 

it' sin (bj 



(1.24) r w i?- i?0tan£ + Qiitan 2 ^ 2 + O (0 3 ) «i?-.R0tan£ 



(1-25) £ 



1 1 i?' 



dr i?tan£ RR-R' 

The expression for t approximates to: 



R R 



The metric elements are 



(1.27) g« = c 2 (l-f) 

(1.28) g rr = -(l-^) 

(1.29) 5ee = 



r - 

„2 



(1.30) <7 00 = -r 2 sin 2 

(1.31) 3*0 = c^asin 2 ^ 

where a is the Schwarzschild radius and a is the specific angular momentum of 
the field source. 

The integral for t is: 

R' R' 
1 R' , aa f 1 , 1 



ft -R 

R-R' a / J?Vi \ , 



Recalling that <^ = fit* 
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(1.34) 
(1.35) 



R-R 1 a ( R'nu \ , R 

l-TTT^^ajlog- 



U 



c c 

R-R 1 a 
h - 

c 



J? (R - R') 

R'flti o\, i? 
1 = I log — 



c c\ R cj R' 
Under conditions where the length of the interferometer arm I -C R 



(1.36) 



1+ R(}- a c 



This expression reproduces the Shapiro time delay || with a small correction 
due to the angular momentum of the field source. 

For the return trip to the center of the interferometer, 



(1.37) 



R> 



cos cf> — tan C sin <f> 

In this case both <f> and t are zero at the starting point i.e., the end mirror. Once 
again 

Rcos6„-R' R-R' 



(1.38) 
(1.39) 



tan ( = 



dr 



i?sin( 
1 

i?'tanC 



R<t> e 



The time for the return trip (using now a positive sign for the radical in (1.18)), 
t e ~ ij. The round-trip takes 



(1.40) 



I 

if — ■ tp 2 

c 



a / „ a 
^R-i 1 -^ 



We note that the corrections add in the round-trip travel time. 
In terms of proper time 



(1.41) 



(1.42) 



-\j 9tt + 2gtt/>Q + g^ffltr 

a aaRVL f RQ\ 2 I 
2- 

c 



1 + 

R R 2 c 



R\ cJ 



We define three dimcnsionless parameters: 



(1.43) 
(1.44) 

(1.45) 







R 

a 

R 
QR 



Keeping terms which are of lowest order in the parameters fj,, n and f3, the radial 
proper time of flight is 
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(1.46) 



2 l - 

c 



There is a general relativistic contribution of ~ a/R. 

Substituting the metric elements in the expression for r e and carrying out 
the calculations to the same order (with I ~ 



(1.47) 
(1.48) 



~ 2 




(hk) + (1 - /z) 
-fj,+ 2fiK/3 - (f 



Note that there is no first order general relativistic contribution. However, in the 
difference in the time of flight between the two arms 



(1.49) St ve = T r — T, 

there is a general relativistic contribution. 



1.3. Light beam traveling along the 6— direction. The light ray travels from 
the beam splitter along the local meridian, either towards the North or South. In 
this case dr ~ 0. The travel time from the beam splitter to the end mirror is 



(1.50) 



9tt dx J gtt \ 
o \ 



9t4>\j^) -<lu\!lnn 



where the angle x is complementary to 0. The integral limits are chosen as and $ 
instead of x^nd x 2 because $ is the angle subtended by the interferometer arms. 
The <7^„'s maintain the values they have at the equator because we approximate 
the rays as straight lines on a sphere of radius R. The rays have a space trajectory 



(1.51) 

The time of flight is 



= kx 



(1.52) 

(1.53) 
Also 



tN = 

gtt 



gu 



9t4, k2 ~ 9tt (gee + 



tN - —V-gttgee$ - — fc* 
gu gu 



(1.54) 
(1.55) 



k 



n 
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Thus 
(1.56) 



t N = VzMg* 



9u + gt<pQ 

The reverse path, calculated with a (— ) sign outside the radical in (1.51), is just 
the same i.e., t$ = ijv, thus the total time of flight is 



(1.57) t M 



t s + t N = 2 V - 9tt9ee ^ 



9tt + 3*0^ 



(1.58) t M ^ 2- 
In terms of proper time 

(1.59) r m 

(1.60) T m = 2- 



7^ 



t M \Jl - h + 2iik[3- [3 2 



1 



Notice again the absence of general relativistic contribution; this holds true for light 
travel time on the tangent or 6 — <fi plane. 
The differences in proper times are 



(1.61) 
(1.62) 



St 

u 1 em 



= T R - Tr 



C 



-J 1 



Expressions (1.61) and (1.62) suggest the experiment. 

It is worth emphasizing that the differences are in proper time; the result is 
independent of any coordinate system. Indeed a careful calculation using isotropic 
or harmonic Schwarzschild coordinates confirms this assertion as does a calculation 
using an expansion about a point. 

These results are not really new; they conform to other relativistic time de- 
lay phenomena such as the radar-echo experiment of Shapiro as well as the clock 
corrections needed for GPS navigation satellites. 

1.4. Numerical values. For the Milky Way Galaxy (assuming it is a homogeneous 
disk): 



(1.63) 


a 
R 


o^Tl -lO-tolO- 
2.8 x 10 20 m 


(1.64) 


[3 ~ 


1CT 3 


(1.65) 


a = 


R2il = R Vp ~ (10 14 to 10 18 ) m 


(1.66) 


a 
R 


1CT 6 to 1(T 2 



for peripheral velocity v p ~ 600km/s. 

For an interferometer located on Earth of length 10cm <5r rm ~ 6 x 10~ 15 seconds. 
This corresponds to an apparent increase in length due to the Galactic center of 
~ 1000 A in the radial arm. 

Corresponding figures for the Sun-Earth system are: 
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(1.67) 
(1.68) 



R 





a 



10 3 ? 
10" 



'm 



rn 




~ 10 



-4 



And at Earth's surface: 



(1.69) 



a 



10" 2 TO 



~ 10" s ; 



(1.70) 



R 

(3 ~ 



6 x 10 6 m 
10~ 6 



The Galaxy is the major source of metric perturbation in the vicinity of Earth. 

Note that the Galactic influence may need to be built into the clock rate adjust- 
ment in the GPS navigation system; currently only the much smaller Earth's effect 
is built-in. 

Why is there a measurable effect proportional to the potential a/r, when one 
can always define a new set of coordinates to make a/r = at r? The answer lies 
in the extended reach of the interferometer; it is not possible, except in flat space, 
to make the potential zero everywhere over an extended region. The interferometer 
measures a potential average over a 10cm x 10cm x 10cm space-like region. This 
is not zero; and cannot be made zero everywhere within the region using a single 
transformation. 

An alternative explanation is this: Two null vectors are transported simultane- 
ously, in closed paths, along two orthogonal axes; when they return to the origin 
there is an angle defect between them, which is, as expected, a first order effect in 
a/r. 

1.5. Experiment, Apparatus and Noise. As the interferometer rotates with 
Earth the output measures alternately, ST rm and Sr em . The signal appears as 
alternating bright and dark regions in the combined beam. The interferometer is 
located at the equator (assumed only to simplify the discussion) , with the two arms 
aligned West-East and North-South. As the interferometer rotates with Earth, the 
North-South arm maintains its relative alignment (this is not strictly true but we 
make this assumption to simplify the discussion) with the Galactic center while the 
orientation of West-East arm alternates between the radial and azimuthal directions 
every 12 hours. The signal is sinusoidal, modulated with a period of 12 hours; it 
can be recovered using a homodyne detector. This is a null experiment in the sense 
that the appropriate phase of the dynamic output is zero in the absence of any 
general relativistic effect, the output measures only deviations from a flat metric. 

It is worth emphasizing that the computed general relativistic fringe shift is for 
each round-trip of the light beam. The fringe shift is the result of transporting, in 
a closed loop, a null vector in curved space. Each subsequent round-trip adds to 
the fringe shift from the previous cycle. By contrast a fringe shift in flat space is 
not cumulative; it is static. 

Each round-trip, which takes 0.67 x 10~ 9 scconds, contributes 1/5 of a fringe 
shift. Thus the fringe shift is wiped out every 5 round-trips or in 3 x 10~ 9 scconds. 
In order to observe the fringe shift one needs to strobe the output signal with a 
frequency which is the inverse of 0.67 x 10~ 9 seconds or 1.5 GHz. 
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Anisotropic Lorentz contraction of the interferometer arms (1-47) will occur with 
a magnitude tx (1/2) (3 w 10~ 7 . This is smaller than the expected signal but also 
90° out of phase with it, hence distinguishable. 

Although we have assumed that the length of both interferometer arms is exactly 
equal, this will not be so in practice, nor is it necessary. Any remnant inequality 
will show up as a zero-offset in the sinusoidal signal. 

A time delay of 10 -16 sees, corresponds to ~ 1/5 of a fringe shift. Is this 
measurable in the presence of noise from sources such as mechanical/thermal noise, 
photodetector noise; mechanical distortions such as sag and tilt, seismic noise, 
laser intensity fluctuations, laser wavelength/frequency fluctuations etc.[Q? The 
effect we are looking for is independent of wavelength, so neither laser wavelength 
nor intensity fluctuations nor drift will affect the result if the apparatus is kept 
in a vacuum environment. However, the interference condition does depend on 
wavelength, so the laser wavelength needs to be stable to <C 1/5 of a wavelength. 

Among the noise sources listed, the most relevant is thermal stability since we are 
looking for a change in apparent length. Low-expansion substances such as Zerodur 
or sapphire would be suitable platform materials. With expansion coefficients of 
~ 10~ 6 a length of 10cm can be kept to within -C ±1000 A by controlling the 
temperature to within ±10 _3 K. This is a modest challenge at room temperatures. 

Vibration isolation is effective with increasing resonance frequencies; this also 
makes sapphire or Zerodur suitable materials because of their mechanical stiffness. 

Integration intervals over several days or weeks will require reference oscillators 
stable over long periods. Oscillators are available with fractional frequency stabili- 
ties of ±10 -16 over months. 

One may ask why this relatively large effect was not detected in experiments 
going back to that of Michelson and Morley. We have studied published accounts 
(2) and have concluded, tentatively, that the apparatuses and methods used were 
neither sensitive to nor designed to measure what we are proposing. 

1.6. Proposed experiments. Two experiments are proposed, and a recommen- 
dation: 

1) Measure mass of Galaxy ( the part within the Galactic center and Solar 
system) 

2) Possible improvement in the uncertainties in the Shapiro time delay using the 
effect of the Sun on the interferometer. 

3) An immediate recommendation: Make clock rate corrections to satellite based 
navigation systems. 

The second would be a challenge; it may require a cryogenic environment where 
temperature stabilities of ±10 _6 K are possible. This, combined with an expansion 
coefficient of sapphire at 4K which is estimated to be ~ lO -11 -^ -1 , may reduce 
thermal stability to ~ 10~ 17 |). 

The Galactic center would provide a strong background for the second experi- 
ment. It may be possible to extract the solar signal by monitoring the interferometer 
output during months of June and December when the Sun, Earth and Galactic 
center are aligned, then during March and September when they form a right an- 
gle triangle. The two signals, with different amplitudes, will have a relative phase 
which will progress from 0° to 90° between December and March. 

It would appear that the Galaxy can be weighed with a table-top device which 
is sensitive to metric perturbations over cosmic distances. 
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